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ON FIGURES OF EQUILIBRIUM OF A ROTATING COMPRESSIBLE 
FLUID MASS; CERTAIN NEGATIVE RESULTS* 

BY 

E. J. MOULTON 

1. Introduction. The problem of determining from theoretical considera- 
tions the shape of celestial bodies has led to extensive discussions of figures of 
equilibrium of a rotating fluid mass. Especial difficulties have been en- 
countered when the fluid is assumed compressible, the case in nature. They 
are in fact so great that various writers have considered it worth while to 
establish that certain figures, for a certain type of law of density of the fluid 
mass, are not figures of equilibrium. Thus Hamyf proved that an ellipsoidal 
mass composed of a finite number of layers of homogeneous fluid in which 
the density increases layer by layer from surface to center and the bounding 
surfaces of the layers are ellipsoids coaxial with the surface ellipsoid, cannot 
be a figure of equilibrium. VeronnetJ extended Hamy's theorem to the case 
in which the density is any continuous function increasing from surface to 
center. Volterra§ arrived at the same negative conclusion, supposing the 
equidensity surfaces to be homothetic ellipsoids, and the density function to 
be any positive, finite, integrable function not a constant. 

In this paper other negative theorems are established. Surfaces of great 
generality are considered, as contrasted with the ellipsoidal surfaces discussed 
by earlier writers, but generality in the conclusion is lost in so far as there 
may be for a given surface a finite number of values of a certain compressi- 
bility parameter for which the conclusion is not true. The surfaces discussed 
include all ellipsoids, and the results for ellipsoids are new since it is not 
assumed that equidensity surfaces are ellipsoids. 

2. Notation and fundamental equations. A fluid mass is considered whose 
surface is S and density at any point P is p; such a mass is referred to briefly 

* Presented to the Society, under a slightly different title, December 27, 1913. 

t Thesis (Paris), Etude de la figure des corps celestes, 1887. Reproduced by Tisserand, 
Traile de Mecanique Celeste, vol. 2, pp. 186-192. 

I Thesis (Paris), Rotation de I'elKpsoide heterogene et figure exacte de la Terre, 1912, Journal 
de Mathematiques, ser. 6, vol. 8 (1912), pp. 331-465. 

§ Sur la stratification d'une masse fluide en iquilibre, Acta Mathematica, 27 (1903), 
pp. 105-124. 
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as the mass ( S , p ) . The mass is rotating with an angular velocity w about 
an axis A; the angular velocity is constant throughout the mass. The only 
forces acting arise from the mutual attractions of the particles under the 
newtonian law, and from the internal pressure. 

In order for the fluid mass to be in equilibrium the following conditions 
must be satisfied: 

I. At every point of the fluid 

(1) ■ dp=fpdU, 

where dp and dU are the differentials of the pressure and of the potential of 
the attraction and centrifugal force, and where / is the gravitation constant. 

II. On the surface S the potential U and the density p are constant. 

III. The center of gravity of the mass ( S , p ) lies on the axis A of rotation. 
If, for a mass ( S , p ) rotating with an angular velocity w about an axis A , 

conditions I, II, III are satisfied, then ( S , p , co , A) will be called a system 
of equilibrium. 
Turning to the equations of the problem, we have first 



(2) ^ = ^+J 



2/ ' J A 

where R is the distance from the axis of rotation A to the point P at which 
U is computed, and where A is the distance from P to the volume element dr , 
and the integral is extended throughout the mass ( S , p ) . Next we assume 
a relation between the pressure and density in the mass to hold, expressed 
by the equation 

(3) hdp = Jp ~ dp , 

where A is a compressibility parameter, and F{p, h) is a polynomial* in p, 
is analytic in h in the vicinity oi h = 0, and is such that 

op Jp=p„ 

if Po > , po being a constant. 

Equations (1) and (3) are satisfied if and only if 

dF(p, h ) J 

;, dp = hdU or p = v . 

op 

The latter equation may be discarded in our problem. The solutions of the 
former under the condition that p = po > at a fixed point Po in the mass 
( S, p) are considered; they are solutions of 

(4^ F{p^,h) =F{po,h)+h{U^-Uo), 

* The function F (p , h) may be greatly generalized in the later argument; see § 7. 
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where subscripts are used to indicate the point at which the functional values 
are taken, Pi being a variable point in the mass ( S , p ) . In view of equation 
(2), we may write (4) in the form 

7 2 /* 

(5) F{puh) =F{po,h)+^{Rl-Rl) + hj K.^p.dr,, 

1 1 

where A12 is the distance from Pi to the volume element dT2, and A02 has a 
similar meaning. 

Equation (5) is the fundamental integral equation which the density function 
must satisfy. It has its simplest form when F (p, h) = p. In this case (3) 
is the Laplace pressure-density equation* 

(3') hdp = fpdp; 

and (5) may be written 

(5') p'l = 1 + hqim - Rl) + hfKup'2dr2, 

where 

(5' a) p'l = pi/po, q = wV2/po. 

3. More exact statement of the problem. The general problem of deter- 
mining figures of equilibrium from our point of view may now be stated thus : 
Given an axis of rotation A and a pressure-density equation of the form (3), 
it is required to find an angular velocity o) , a surface S , and a constant po > , 
such that a solution p of equation (5) makes ( S , p , co , ^ ) a system of equi- 
librium. Our method will be to make the compressibility parameter h play 
a fundamental role; that is, various functions will be expanded as power 
series in this parameter. In this paperf we shall let S be any one of an ex- 
tensive class of surfaces and Po any interior point, and consider the possibility 
of determining w and po as functions of h such that when equation (5) is satis- 
fied CS , p , 0}, A) is a system of equilibrium for some continuous range of 
values of h . The negative conclusion we shall reach shows that to solve the 
problem by our method it is necessary to let S vary with h . 

In order to describe S accurately, we make the following definitions. By a 
regular surface we understand one which can be cut into a finite number of 
pieces, each of which has at every point a continuously turning tangent plane; 
by a simple surface, one which is continuous and has no multiple points or 
lines. The surface S ivill be chosen from the class of closed, simple, regular surfaces. 

We shall assume that co and po can be determined as desired and derive 

* The pressure-density equation dp = hpdp was used by Lapljice (Mecanique Celeste, 
book XI, chap. IV); it is to be observed that our h is proportional to the reciprocal of the 
Laplace h . For us, h = gives an incompressible fluid. 

t In a later paper I propose to adopt a different method, and let S vary with h . 
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necessary consequences. We first consider the special case when (5) reduces 
to the simpler equation (5'). 

4. The solution of the integral equation (5'). This is a linear integral 
equation of the Fredholm type, but presents an unusual difficulty in the fact 
that the kernel ku becomes infinite not only when the variable points Pi 
and P2 coincide, but also when P2 coincides with the fixed point Po ■ This 
difficulty is met by considering the iterated equation* 

(7) Pl = m + h^ fl^l2p2dT2, 

where we have dropped accents and written 

,71 = 1 + hqim - Rl) + /^ J«i2 [1 + hq{Rl - Rl)]dr2, 

(8) ^"2 = I Ki3 K32 rfrs 

~J LA13A32 ~ A03A32J ^^~Ao2./ LAi3~Ao3J ^*' 

The function rn is readily seen to be continuous in Pi throughout S; more- 
over the last two integrals are continuous functionsf in Pi and P2 throughout 
S , and hence we may write { 

(2) I ' 

K12 = M12 + M12 • T~ , 
ii02 

where the only discontinuity comes in through the factor I/A02 • If we use 
polar coordinates with the pole at Po , we have 

A02 = ^2 , dT2 = rl sin 62 dr2 dB2 d<j)2 , 

and hence the function 771 and the kernel i^ii rl sin 62 are continuous in 
{ri,di,<j)i) and {ri, 61, (f>i, r2, 62, (f>2) respectively. Thus the apparent diffi- 
culty in equation (7), due to the way in which the kernel becomes infinite, 
disappears. 

We know then from the general theory of integral equations,! that equa- 

* That is, the equation obtained from (5') by changing the subscripts 1 and 2 to 2 and 3 
and substituting the right member for p2 in the integral of (5'), and then altering the notation 
slightly. It is here, where an interchange of order of integration occurs, and later in this 
article in connection with solving (7), that some assumption concerning the character of S is 
chiefly needed. It could be of less restricted class. 

t The continuity of the function given by the first integral follows from the inequalities 
given by Fredholm, Acta Mathematica, vol. 27 (1903), p. 387. The second inte- 
gral is the potential for a constant density, which is known to be continuous. 

I As a matter of fact /x'u is independent of P2. 

§ See, for example, B6cher, An Introduction to the Study of Integral Equations; or Lalesco, 
Introduction a la Theorie des Equations Integrates. 
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tion (7) has, except for isolated finite values oi h, a, unique solution continuous 
in Pi , and it readily follows that the same is true of (5'). Also we may write 
the solution 

(9) Pl = Vl + K'fK?i7l2dT2, 

where Kfj^ is, except for isolated finite values of h, analytic in h and con- 
tinuous in Pi and P2 throughout S . In particular, Kf^ is analytic in h at 
^ = 0. 

By use of equation (8) we write 

(10) pi = 1 + hqim - Rl) + hSi{l) + ¥ qSiiR' - Rl) , 
where, if / is any function whose value at P» is fi , 

Slif) = j'KnJ2dr2 - /j/Kfi[/2 + }iSK..zfzdrz'\dT2. 

We observe that when/ is independent of h, Si{f) is analytic in h except 
for isolated finite values of h , and in particular is analytic a,t h = . 

5. The first conclusion. .If co and po can be determined so that when pi 
is a solution of (5') (S,p,w,^)isa system of equilibrium for other than 
isolated finite values of h , then it is possible to determine q{ = w^/2/po ) so 
that the function pi given by (10) is effective in satisfying equilibrium con- 
ditions II and III. In the first place q may be so chosen that the values of p 
at any two points Pi and P2 on S are the same. The value of q for this 
purpose is readily found to be 

(11) ,= SAD -S.il) 



Ri - Rl + h[ Si{R! - Rl) + S2{R', - Rt)]' 

where subscripts on the S's indicate the points at which they are computed. 
Since g ( = w^/2/po ) is constant as far as Pi and P2 are concerned, its charac- 
terization as to h is obtained by choosing Pi and P2 so that i?i #= P2 • We 
thus find that q as given by (11) is analytic in h except for isolated finite values, 
and in particular is analytic at h = 0. This value of q, therefore, when 
substituted in (10), determines a function p which is in general analytic in h, 
in particular in the vicinity oi h = 0, and which for all such values of h is 
constant on S . Also since the center of gravity of ( S , p ) is on the axis A 
for other than isolated finite values of A , it is in general on A , and in particular 
for all values of h near h = 0. 

The argument thus far reduces the consideration of values of h to those in 
the vicinity o( h = 0. We now write 

(12) Pi = t,p['^hi; q^t^q^^n^, 
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where p^/' is continuous throughout S . Equilibrium conditions I and II 
require (a) that p'-p be constant on S , and (6) that the center of gravity of {1, , p'^^) 
lie on the axis A . 

By substituting (12) into (5') and equating coefficients of Hke power of h, 
we find that 

^^^^ p('' =?('-» (^f-iiS)+ /«12Pr'^^T3 (i=l,2, 3,.--). 

The right member of (13) is, except for an added constant, the potential at 
Pi, for a certain angular velocity, of the mass (S, p^'"^^) . In view of the 
italicized statement in the preceding paragraph, taking i = 1 , we now have 
the following negative result. 

Theorem 1. Let an axis A of rotation and a closed, simple, regular surface* 
S he chosen arbitrarily. Then, {a) if S is not an equipotential surface of the 
mass (2,1) when rotating with a suitable constant angular velocity, or (b) if 
the center of gravity of the homogeneous mass ( S , 1) is not on A , there exist only 
isolated finite values of h for which « and po 'may be so chosen that, if p is a con- 
tinuous solution of (5'), C^ , p , d) , A) is a system of equilibrium. 

6. The case of an ellipsoid. If S is an ellipsoid and A is a diameter, the 
center of gravity of the mass ( S , 1 ) lies on A , but S may or may not be an 
equipotential surface of (S, 1) when rotating with a proper angular velocity; 
if it is not, Theorem 1 applies; if it is, we proceed one step further. 

For this purpose, we use the following method of computing the newtonian 
potential. Let 2 be 

(14) -.+!+-. = 1> 

and let us write 

^2 y2 ^2 yS yi ^2 

^ " -^ ~ a^ ~ 62 - -2; J' = 1 - ^_^- - 52 _^ 5 - J-j^-g • 

Then if the density is constant on the surface p. = const., the newtonian 
potential is given byf 



(15) 






'^(a^'+s){b^ + s){c^ + s) 



irabc 
Now suppose Po is taken on the axis A , so that Po = . Taking i = 1 



* A more extensive class of surfaces might be used; e. g., S might be composed of a finite 
number of non-intersecting closed, simple, regular surfaces, and the conclusion still be true. 

t This formula was used by Volterra, loc. cit., p. 108, with a slightly different jiotation. 
The proof is readily made by comparison with pp. 138-140 of Mathieu's Theorie du Potential, 
where a similar formula is derived for the potential at an exterior point. 
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in (13), we have then from (15) 

(16) p?> = gC' m - if 1 x! - if 2 yi -Mszl, 

where 

ds „ r°° ds 



^' = i ra'^+Trmj)' ^^' = i 



M 



/'OO 



ib' + s)Dis)' 
ds 



i(^ + s)D(s)- 



It is readily proved that since pi^ is constant on Si, and hence the right 
member of (16) is constant in xi and yi when zi is eliminated by use of (14), 
the axis A of rotation must be a principal axis of S . Suppose A is the s-axis. 
On eliminating zi from (16) and equating the coefficients of xl and ?/? to zero, 
we have 

2 2 

gW) - Mi-~M, = 0, </("> - M2 - pifs = 0. 

The ellipsoid must therefore have axes so related that for some value of q'-"^ 
these two equations hold simultaneously.* For such ellipsoids 

(17) pV>= -c^if3[^; + | + ^] = c^if3(M-l). 

We consider next p'^' . From (13) we have 

pf = g«iif+/Kl2P?'^T2. 

The integral is, except for an added constant, the newtonian potential for the 
density function pi". It may therefore be computed by use of (15) and (17). 
It is seen that pf contains terms of second and fourth degrees in xi, yi, zi; 
the latter are 



Jo V a^ + 5 



¥ + s ' c' + sj D{s)' 



When zi is eliminated by use of (14) the fourth degree terms must vanish; 
hence 

p (g^ - (?)sds _ r" (6^ - (^)sds 

Hence we have a = b = c; that is S is a sphere. Hence 

Theorem 2. Even if S is an ellipsoid, not a sphere, and A is a diameter 
of S , the conclusion of Theorem \ is true. 

* These conditions are satisfied if the ellipsoid is of revolution or if it is one of the figures 
of equilibrium discovered by Jacobi — the Jacobian ellipsoids. 
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7. The more general pressure-density equations. In the preceding para- 
graphs our equations come from the initial use of the Laplace pressure-density 
equation (3')- We now consider the consequences of the initial use of the 
more general equation (3). We choose the surface S and axis A as in the 
preceding discussion and attempt the formal determination of po, p , and co^ 
as functions analytic in h in the vicinity of A = such that equation (5') 
and equilibrium conditions II and III are satisfied for all values of h in the 
vicinity of A = . If we write 

(18) PO = E p« h'; p = Z p^'-) h'; ^ = E coC^ h\ 



i=0 



conditions II and III require that p^'^ be constant on S and that the center 
of gravity of the mass (S , p*'^ ) be on ^ . 

If we substitute (18) in (5) and take the value a,t h = oi F (pi, h) and 
its first and second ^-derivatives,* we have 



(19a) F(pf,0)=F(pT,0), 

(196) 



aF(pr.o) ,„ ^ dFipf\o) dF(pT,o) ,,,^dF(pT,o) 



dpi '^' ' dh dpo '^'' ' dh 

+ o,^0HRl-Rl)+fKi2pTdr,, 



dFipT,0) ,„ d^F(pT,0) ,,,. d'F(pr,0) (,) 
d^^-'^+ dpi "' +2 dpidh P' 

d'FipT,0) _ dF{p':\0) ,,, d'Fip<:\0) ^,^, 
(19c) dh' dpo "» + dpi ■"» 

d'F(pf\0) d'F(,^\0) 

+ ^ dpodh "» + dh' +2«^Hiii-i?„) 

+ / Kl2 P*2 ' dT2 . 

Equation (19a) is satisfied by p^"* = p'-g\ which makes pi equal po, as it 
should, when the point Pi coincides with P. Since p^J*^ is a constant, Pi"' is 
constant on S, satisfying one equilibrium condition. The center of gravity 
condition requires that the center of gravity of the mass ( S , p'J" ) , and hence 
of (S, 1), lie on the axis A. Solving equation (196) for p'/^ we have a 
constant added to a constant multiple of 



,(0) 



Rl+fKnpfdT2 



* It is only in computing these derivatives that we need restrictions on the kind of func- 
tion F (f),h) may be. Any function which leads to equations (19a), (196), (19c) on 4ifferen- 
tiating asioh, and for which the coefficient of p^'' and p^^^ is not zero, may be used. 
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This expression is therefore constant on 2 under precisely the same conditions 
as the function p^^ of the preceding paragraphs, and differs from that function 
only by an added and a multiplied constant. Hence when S is an ellipsoid 
and A is a. diameter of S , the function p''^ is either not constant on S or it is 
given by a formula similar to (17) : 

When we consider pj^', as found from (19c), we see that, since p*/' is constant 
on S , p'-2^ is constant on 2 if and only if 

is constant on S . This leads at once, as before, to the conclusion a = b = c. 

Hence we have 

Theorem 3. Let S and A be as described in Theorem 1. Under the same 
hypotheses (a), (b), or the hypothesis of Theorem 2, it is impossible to determine 
CO , po , and p as power series in h so that equation (5) shall be satisfied, p shall 
be constant when h = 0, and ( S , p , co , A) shall be a system of equilibrium for 
all values of h near h = 0. . 



